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AN  IMPROVED  ANALYTICAL  MODEL  OF  THE  INTERACTION  OF  DOMES 
AND  TRANSDUCERS  DURING  TRANSMISSION 

I.  INTRODUCTION 

The  purpose  of  this  technical  memorandum  is  to  describe 
an  improvement  in  the  specification  of  transducer  velocity 
distributions  which  are  used  in  the  analysis  of  dome -transducer 
interactions.  The  mathematical  model  used  in  these  studies 
consists  of  a  circular  cylinder  and  concentric  shell.  Typically 
the  interaction  problem  has  been  studied  by  obtaining  a  solution 
to  the  boundary  value  problem  defined  by  the  geometry  of  the 
cylinder  and  concentric  shell. 

In  previous  studies  a  velocity  equivalent  to  the  radial 
component  of  a  plane  sound  pressure  wave  has  been  specified  on 
a  portion  of  the  cylinder  face,  the  remainder  of  the  cylinder 
being  rigid.  The  portion  of  the  cylinder  face  having  nonzero 
velocity  represents  the  active  elements  of  the  transducer.  For 
such  a  velocity  specification,  the  velocity  (amplitude  and  phase) 
is  continuous  across  the  active  face.  This  specification  is 
somewhat  fictitious  since  it  is  equivalent  to  an  infinite  number 
of  point  sources.  In  a  real  transducer  there  are,  of  course,  a 
finite  number  of  active  elements  having  finite  dimensions  and 
discrete  velocity  amplitudes  and  phases. 

Recently,  a  computational  procedure  has  been  developed 
which  permits  the  active  portion  of  the  cylinder  face  to  be 
separated  into  finite  segments,  each  of  which  has  a  discrete 
velocity  amplitude  and  phase.  The  mathematical  form  of  the 
function  describing  the  velocity  boundary  condition  specified  on 
the  cylinder  has  a  finite  number  of  finite  discontinuities  as 
compared  to  the  previous  continuous  description.  This  new  boundary 
specification  is  clearly  a  more  realistic  one. 
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One  important  result  of  the  improved  boundary  specification 
is  that  dome- transducer  interactions  now  can  be  analyzed  for  various 
transducer  element  shading  and  phasing  factors.  This  generality, 
not  available  in  previous  boundary  specifications,  will  permit 
the  investigation  of  an  "optimum**  element  shading  and  phasing 
including  the  effects  of  the  dome. 

In  addition  to  the  mathematical^de tails  of  the  new  boundary 
specification,  some  computed  results  for  the  near-  and  farfields 
of  the  dome- transducer  are  included  in  this  memorandum.  These 
results  are  compared  with  results  previously  obtained  for  a 
continuous  velocity  distribution.  Results  are  also  presented 
for  several  types  of  shading  to  indicate  dome  effects  on  the 
nearfield  sound  pressure  distribution  and  on  beamwidth  and  side 
lobe  suppression  as  a  function  of  transducer  shading. 
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II.  FORMULATION 


The  interaction  of  the  dome  and  transducer  is  analyzed 
by  obtaining  a  solution  to  the  scalar  wave  equation  for  boundary 
conditions  prescribed  on  the  transducer  and  dome  surfaces.  The 
geometry  of  the  dome- transducer  model  considered  in  this  memo¬ 
randum  is  identical  to  the  geometry  considered  in  previous  studies 
(see  Ref.  1)  and  is  shov/n  in  Fig.  1. 


The  boundary  conditions  typically  employed  in  the  boundary 
value  problem  are  that  the  normal  particle  velocity  of  the  fluid 
have  a  prescribed  value  at  the  transducer  face  and  that  the  normal 
particle  velocity  be  continuous  at  the  shell- fluid  interfaces. 

In  previous  studies  the  normal  particle  velocity  specified  on 
the  active  portion  of  the  transducer  face  was  equal  to  the  radial 
component  of  the  particle  velocity  of  a  plane  pressure  wave,  the 
remainder  of  the  cylinder  being  rigid.  This  boundary  condition 
has  the  mathematical  representation  * 


r=a 


=  cosB  e 


ika  cos 


0  <  i9  I  <  cp 


=  0 


cp  <  1 0  !  '*'  TT 


(1) 


where 

V  is  the  particle  velocity  of  the  fluid  normal  to  the 
transducer 

2tt 

k  =  and  \  is  the  wavelength  of  the  transmitted  wave, 

and  a  is  the  radius  of  the  transducer. 

*The  time  component  e  ,  where  uj  is  the  angular  frequency,  will 
be  suppressed. 


3 


REGION  2 


REGION  I 


DOME 


ACTIVE  STRIP 


TRANSDUCER 


FIG.  I  -  GEOMETRY  FOR  A  CYLINDRICAL  TRANSDUCER 
RADIATING  IN  A  CONCENTRIC  DOME 
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This  boundary  representation  has  the  effect  of  shading  each 
point  on  the  active  portion  of  the  boundary  by  a  cos  0  factor 
and  phasing  each  point  to  a  straight  line.  In  this  memorandum 
a  velocity  condition  is  employed  having  the  representation 

vlr^  =  CVj]  ,  (2) 

where 

[Vj]  is  the  set  of  head  velocities  of  the  active  transducer 
elements,  and  is  the  complex  normal  velocity  of  the  jth  element 
in  the  transducer. 

It  is  assumed  that  the  new  velocity  specification  is  an  even 
function  of  9  although  other  cases  can  be  handled  with  little 
increase  in  difficulty.  A  comparison  of  the  continuous  and 
discrete  boundary  conditions  is  shown  graphically  in  Fig.  2 

The  velocity  of  any  particular  active  element  in  the 
transducer  array  can  be  written  as 

iq) . 

Vj  =  IVjle  J  (3) 

where  jV.|  is  the  velocity  amplitude  and  cp  .  is  the  phase  of  the 
1  J 

jth  element. 

The  velocity  specification  for  a  single  element  of  width 
2P  and  angular  position  Y=(2j-l)a  as  shown  in  Fig.  3  is 

vL_  =0  0  <  |e|  <  Y  -  @ 

icp  . 

=  jv^le  J  Y-P<ie!  <  Y  +  P  (4) 

Y+§< 1 0  I  <  n 


0 


NORMALIZED  ELEMENT  VELOCITY  AMPLITUDE  IVjl 


FIG.  2  -CONTINUOUS  AND  DISCRETE  BOUNDARY  CONDITIONS 
PRESCRIBED  FOR  THE  DOME-TRANSDUCER  MODEL 
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FIG.  3  -  TRANSDUCER  ELEMENT  GEOMETRY 
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This  distribution  has  a  Fourier  Series  representation 


vi  =  y  S  cos  n0 
'  r=a  ZLi  n 

n^O 


(5) 


where 


Y4^ 


(>„  =  rr-  r  '  "  |v|e^  cos  nede 

n  ^  '  Y-P 


=  ^ 


rr 


V\e^^  ,  n=0 


~  Ivle^'^  [sin  n(Y-{-P)  -  sin  n(Y-P)3, 


n>0  , 


and  s  =1  for  n  =  0,  e  =2  for  n>0. 
n  ’  n 


Since  the  problem  is  linear,  the  velocity  distribution 
for  several  elements  can  be  obtained  by  vector  addition  of  the 
velocities  of  individual  elements.  Thus  if  an  array  of  2N 
equally-spaced,  active  elements  is  chosen  as  shown  in  Fig.  3, 
with  the  jth  element  having  a  velocity  given  by  Eq.  3, 
the  boundary  condition  takes  the  form 


ICO. 


''Ir=a  =  V 


=  V^e 


L'Og 


.  ICO  . 

V  e  ^ 
J 


ICO 


V 


0 


N 


a-p<|0|<a+p 

3a-p<|0|<3a-f^ 

(2j-l)a-p<|0|<(2j-l)a+^ 

( 2N- 1 )a-p< 1 0  I < ( 2N- 1 )a+0 
elsewhere 


(6) 


where  2a  is  the  element  spacing  and  2P  is  the  element  face  width. 
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In  a  manner  analogous  to  Eq.  5,  Eq.  6  can  be  given 
a  Fourier  series  representation 


V 


r=a 


I 


n*0 


6  cos  n9 
n 


where 


and 


{sin  nL(2j-l)a-fp] 


(7) 


-  sin  n[( 2j-l)a-p] }  ,  n>0. 

The  boundary  condition  is  now  in  a  form  suitable  for  use  in 
determining  the  sound  pressure  field  inside  and  outside  the  dome. 

The  mathematical  details  of  solving  the  boundary  value 
problem  for  a  continuous  boundary  condition  are  given  in  Appendix 
A.  The  solution  of  a  boundary  value  problem  having  a  boundary 
condition  given  by  Eq.  7  is  obtained  in  exactly  the  same  manner 
by  substituting  Eq.  7  for  Eq.  9>  Appendix  'A,  and  carrying  out 
the  indicated  operations  in  Appendix  A. 
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III.  NUMERICAL  RESULTS 

Numerical  results  were  computed  for  the  sound  pressure  in 
the  near-  and  far fields  of  the  transducer  in  order  to  determine 
the  effect  of  the  dome  on  transmitted  beams.  Results  were  obtained 
for  several  sets  of  amplitude  shading  factors  jVj)  as  given  in 
Eq.  6.  For  all  computed  results,  the  phasing  factors  »j,  given  in 
Eq.  6,  were  chosen  so  that  the  elements  were  phased  to  a  straight 
line.  The  results  were  computed  for  parameters  appropriate  to  the 
AN/SQS-26.  Results  with  and  without  the  dome  are  presented.  In 
each  case  the  results  are  normalized  to  the  maximum  sound  pressure 
in  the  distribution. 

The  near-  and  farfield  pressure  distributions  for  a  trans- 

* 

ducer  having  unity  shading  factors  are  given  in  Figs.  4  and  5. 

The  near-  and  farfield  pressure  sound  pressures  presented  in  Figs.  6 
and  7  are  for  amplitude  shading  factors  having  a  cosine  distribu¬ 
tion.  For  this  case  the  amplitude  factor  |V^|  is  obtained  as 

jVjj  “  cos  (2j  *1)0  .  (8) 

A  comparison  of  the  near-  and  farflelds  obtained  with  continuous 
and  discrete  transducer  velocity  distributions  is  given  in  the 
linear  plots  of  Figs.  8,  9,  10,  and  11.  The  former  results, which 
have  continuous  phasing  and  shading  (cosine) ,  were  obtained  from 
Ref.  1. 


For  the  results  presented,  the  transducer  radius  is  a, 
the  dome  radius  is  b,  the  dome  thickness  is  h  (h"0  indicates  no 
dome  is  present), and  r  is  the  radial  distance  from  the  center  of 
the  transducer  at  which  the  sound  pressure  field  was  computed. 
All  dimensions  are  in  inches. 
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-30"  -20*  -!0*  0*  10*  20*  30* 


FIG.  4 -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  =  96  b=120  r:9€ 

SHADING:  UNITY 


-30*  -20*  -10*  0*  10*  20*  30* 


FIG.  5  -  SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  =  96  b=  120  r  =  FARFIELD 

SHADING:  UNITY 


FIG.  6  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  =  96  b=l20  r  =  96 

SHADING  =  COSINE 


TRACOR,  INC  _ I  DWG  778-  57 

AUSTIN.  TEXAS  X/2S/6S 


RM/L<H 


-30*  -20*  -10*  0*  10*  20*  30* 
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FIG.  7  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0=96  b=l20  r  =  FARFlELD 

SHADING*  COSINE 


FIG.  8  'SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  =  96  b=l20  r2  98 

SHADING  =  COSINE 
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0  20  4  0  60  80  100 
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FIG.  9  “SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  =  96  b  =  l20  r=  96 

SHADING  =  COSINE 
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0  10  20  30  40  50 

RELATIVE  BEARING -DEGREES 


FIG.  10  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE:  A  CONCENTRIC  DOME 

0  =  96  b*l20  r  =  TARFIELD 

SHADING  =  COSINE 


TRACOR,  IHC _ I  DWG. 773 -51 

AUSTIN,  TEXAS  *  rt\ aH/LtK 
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Results  for  amplitude  shading  factors  which  have  a  cosine- 
squared  distribution,  i.e., 

|v.|  =  cos^(2j-l)a  ,  (9) 

are  given  in  Figs.  12  and  13,  Shading  factors  corresponding  to 
those  for  the  AN/SQS-26,  Model  XM-2,  were  used  to  compute  the  re¬ 
sults  shown  in  Figs.  14  and  15.  Finally,  shading  factors  derived 
by  the  method  of  Dolph-Tchebyscheff  were  used  in  obtaining  the 
results  shown  in  Figs.  16  and  17. 

A  compilation  of  the  shading  factors  is  given  in  Table  1. 
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FiG.  12  “fOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 


0  =  96 

SHADING  =  COSINE' 


bilEO 


r=  96 


TRACOR,  INC _ I  DWG.  778- 53 
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-30*  -20*  -lO*  0*  10*  20*  30* 


FIG.  13  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

3=96  !J=I20  r=FARFIELD 

SHADING  =  COSINE^ 
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AUSTIN,  TEXAS  8/23/8S  IM/LCH 


FIG.  14  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0=96  b:  120  r  =  96 

SHADING:  AN/SQS -26  ,  XN-2 
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-30*  -20*  -10*  0*  10*  20“  30* 


FIG.  !5  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0=96  b=l20  r=FARFIELD 

SHADING  =  AN/SQS  26  ,  XN2 
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FIG.  16  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

0  :  36  120  r  =  96 

SHADING  =  DOLPH-TCHEBYSCHEPr 


24 


TRACOR,  iMC _ r  DWG.  778-54 

AUSTIN,  TEXAS  S/tS/SB  8M /UH 


-30*  -20*  -10*  0'  10*  20*  30* 


FIG.  17 -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME 

a  =96  b=l20  r:  FARFiaO 

SHADING  =  DOLPH-TCHEBYSCHEFF 
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ELEMENTS 


I 


SHAD  I] 

NG  TYPE  1 

ELEMENT  NQ 

UNITY 

COSINE 

COS^ 

XN-2 

DOLPH 

1 

1.00 

.999 

.999 

1.00 

1.00 

2 

1.00 

.991 

.982 

1.00 

.980 

3 

1.00 

.976 

.953 

1.00 

.925 

4 

1.00 

.954 

.910 

1.00 

.850 

5 

1.00 

.924 

.854 

1.00 

.760 

6 

1.00 

.887 

.737 

1.00 

.665 

V 

1.00 

.843 

.711 

1.00 

.570 

8 

1.00 

.793 

.629 

.933 

.475 

9 

1.00 

.737 

.543 

.708 

.390 

10 

1.00 

.676 

.457 

.631 

.315 

11 

1.00 

.609 

.371 

.562 

.265 

12 

1.00 

.537 

.288 

.501 

.315 

TABLE  I -ELEMENT  AMPLITUDE  SHADING  FACTORS 
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IV.  CONCLUSIONS 

A  computational  procedure  has  been  developed  which  permits 
the  specification  of  a  discrete  velocity  amplitude  and  phase  for 
finite-width  transducer  elements  in  the  dome- transducer  model. 
Previously  it  was  necessary  to  assume  a  continuous  velocity 
distribution  at  the  transducer  in  analyzing  dome- transducer  inter¬ 
actions  . 


A  cor.parison  of  computed  near-  and  farfield  sound  pressures 
indicates  that  the  newer,  more  realistic  transducer  boundary 
specification  results  in  some  changes  when  compared  to  previous 
computations  with  a  continuous  velocity  specification.  For  cases 
with  and  without  the  dome,  the  newer  specification  results  in  a 
rather  high  pressure  amplitude  on  the  inactive  (rigid)  portion  of 
the  transducer.  This  result  indicates  that  there  is  some  inter¬ 
action  among  the  active  transducer  elements  and  the  adjacent  in¬ 
active  elements.  This  effect  is  shown  in  Figs.  8  and  9.  In  the 
farfield  the  side  lobes  are  somewhat  lower  for  the  discontinuous 
boundary  with  and  without  the  dome.  (See  Figs.  10  and  11  .) 

The  purpose  of  shading  a  transducer  array  is,  of  course, 
to  improve  the  directivity  of  the  array.  Other  factors,  however, 
should  be  considered  in  choosing  an  array  shading.  Among  these 
are  smoothness  of  the  nearfield  pressure  and  source  level  - 
decreasing  the  velocity  amplitudes  of  some  of  the  elements  will 
reduce  the  power  transmitted  to  the  water.  The  first  factor  has 
implications  in  regard  to  cavitation  and  element  loading,  each  of 
which  can  affect  system  performance.  The  results  in  Figs.  4-17 
illustrate  the  near-  and  farfield  effects  of  the  dome  for  several 
types  of  element  shading. 

In  regard  to  smoothness  of  the  nearfield  there  does  not 
appear  to  be  much  difference  in  the  results  for  the  shading  factors 
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used  in  the  computations;  the  nearfield  sound  pressure  distribu- 
tions  for  each  type  of  shading  are  all  nonsmooth  with  the  dome  in 
place.  (See  Figs.  4,  6,  12,  14,  and  16.) 

The  results  presented  in  Figs.  5,  7,  13,  15,  and  17  indicate 
how  shading  can  be  employed  to  improve  the  directivity  of  the 
transducer  and  increase  side  lobe  suppression  in  the  farfield. 

An  interesting  result  shown  in  these  figures  is  that  in  cases  of 
extreme  shading,  for  instance  Dolph-Tchebyscheff  shading  shown  in 
Fie  17,  the  increase  in  side  lobe  level  due  to  the  dome  is  more 
than  for  cases  having  lesser  amounts  of  shading.  With  Dolph- 
Tchebyscheff  shading,  the  dome  results  in  an  increase  in  side 
lobe  level  of  approximately  10  dB,  while  for  cosine  shading 
(Fig.  5)  the  side  lobes  are  increased  only  about  5  dB.  Hence,  it 
appears  that  dome  effects  are  more  pronounced  as  the  degree  of 
shading  is  increased.  It  can  be  concluded  that  any  program  aimed 
at  "optimizing*'  beam  patterns  by  shading  the  transducer  elements 
should  include  a  consideration  of  the  effects  of  the  dome. 

It  is  not  possible  at  this  time  to  derive  in  a  closed 
mathematical  form  an  optimum  set  of  element  phase  and  amplitude 
factors  which  include  a  consideration  of  the  dome.  However,  the 
techniques  described  in  this  memorandum  can  be  used  in  developing 
"optimum"  element  shading  including  dome  effects  through  an  inter¬ 
active  procedure. 
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APPENDIX  A.  THE  SOUND  PRESSURE  FIELD  PRODUCED  BY  A  TRANSDUCER 

RADIATING  IN  A  CONCENTRIC  DO^£E 


The  mathematical  model  illustrated  in  Figure  1  consists 
of  an  inf ini te “length  cylinder  and  concentric  shell.  A  velocity 
distribution  independent  of  the  axial  direction  is  specified  on 
the  surface  of  the  cylinder,  and  it  is  desired  to  find  the 
resulting  sound  pressure  field  in  the  fluid  in  Regions  1  and  2 
of  Figure  1.  The  boundary  value  problem  as  stated  is  a  three 
region  problem  (two  fluid  regions  and  the  shell).  In  general 
partial  differential  equations  must  be  solved  both  in  the  two 
fluid  regions  and  in  the  region  of  the  shell  material,  and 
boundary  conditions  must  be  matched  at  shell-fluid  interfaces 
as  well  as  at  the  cylinder-fluid  interface.  However,  the 
problem,  is  considerably  simplified  by  using  the  results  of 
elastic  shell  theory  to  describe  the  motion  of  the  shell's  middle 
surface,  thus  eliminating  the  necessity  of  solving  the  full 
partial  differential  equations  of  linear  elasticity  for  the  shell 
material.  Further,  since  the  shell  is  thin  relative  to  other 
geometrical  parameters, the  boundary  conditions  at  shell-fluid 
interfaces  can  be  matched  at  the  middle  surface  of  the  shell. 

A.  Sound  Pressure  Field  in  the  Fluid 

It  is  assumed  that  in  Regions  1  and  2  the  acoustic 

3 

pressure  p  obeys  the  scalar  wave  equation 


r  Br 
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^  hr' 
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where  r  and  6  are  the  space  coordinates}  t  is  the  time,  and  c  is 
the  propagation  velocity  in  the  fluid. 

For  a  harmonic  time  dependence  the  general  solution 
to  Eq.  1  is^ 
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P  =  y  [A„  H^^^(kr)  +  B„  H^^^(kr)]  cos  nO  (2) 

1111  1111 

n=»0 

where  and  are  the  nth  order  Hankel  functions  of  the 

n  n 

first  and  second  kind,  respectively,  uu  is  the  angular  frequency, 
k  =  uu/c,  and  and  B  are  constants  which  will  be  determined 
by  boundary  conditions.  Equation  2  involves  only  cosine  terms 
as  it  is  anticipated  that  problems  to  be  analyzed  will  possess 
symmetry  about  6=0.  In  Region  1  the  sound  pressure  field  is* 


Pl  =  ^  [A^  H^^Ckr)  +  H^^^(kr)J  cos  nG  .  (3) 

n=0 


The  sound  pressure  and  the  radial  component  v^  of  the  fluid  velocity, 
with  harmonic  time  dependence,  are  related  by 


p  ar  dt 


-iujv^ 


At  the  shell -fluid  interfaces  this  component  of  the  fluid  velocity 
must  equal  the  velocity  of  the  shell.  (This  restriction  does 
not  apply  to  velocity  components  tangential  to  the  interface 
because  the  fluid  is  regarded  as  inviscid  and  hence  will  not 
support  shear  forces.)  Denoting  the  velocity  on  the  cylinder 
by  and  the  radial  component  of  the  velocity  of  the  shell  by  w, 
the  interface  conditions  at  r=a  and  r=b  are 


3Pl 

iujp 

&r 

1 

api 

iu)p 

r=a 


r=b 


=  V 


=  w 


(5) 


(6) 


*The  time  component  e“i-i«t:  will  be  suppressed  in  the  following. 
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In  the  case  of  interest,  a  section  of  the  cylinder  of 
angular  extent  2cp  is  pha>'ed  to  represent  a  plane  velocity  dis¬ 
tribution  of  uniform  amplitude,  while  the  remainder  of  the 
cylinder  is  rigid.  This  condition  can  be  written  as 
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~  cos 
pc 


Q  ^ikacos  8 
6  e 


1®I«P  , 


=  0 


I  0l>5p  . 


(7) 


^  )c]r  c  o  8  @ 

A  Fourier  series  for  the  function  cos  0  e  may  be  obtained 

4 

by  differentiating  the  identity. 


Ikrcos  8  ,  y  £  1"  j  (kr)  cos  ne 

^  n  n 

n=0 


where  e  =1,  c  =2  for  n  >  o,  and  is  the  Bessel  function.  Then 
o  n  n 


n6 


I  0|<cp 


n=0 


=  0 


1  0l>tp  , 


(8) 


where  the  prime  indicates  a  derivative  >7ith  respect  to  the 
argument.  Equation  8  can  be  represented  by  the  Fourier  series 


V  =  ^  6  cos  n0 

r  ipc  n 

n=0 


where 


6_  =  e  /  ^tn 

n  n  y  I  in  m 
m=0 


(9) 


and 
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The  substlcucion  of  Equations  2  and  9  into  Eq.  5  results  in 

7^  /  [A„  H^^^’Cka)  +  '  (ka)]cos  nG  =  7^  /  6  cos  nG .  (10) 

lujp  i—t  *-00  nn  ''j  ipc'  £—*  n 

n=0  n*0 

Noting  that  the  cosine  functions  are  linearly  independent  and 
that  k/oj  =  l/c,  Eq.  10  implies 

H^^^’(ka)  +  H^^^'(ka)  =  .  (11) 

The  radial  component  of  the  shell  velocity  can  be  represented  in 
the  form  ^ 

«D 

w  =  ^  cos  n0  .  (12) 

n=0 

The  boundary  condition  given  by  Eq .  6  then  becomes 

<B  m 

"tV"  X  ^n  '  (kb)]cos  nG  =  ^  cos  nG  (13) 

n=0  n=0 

or 

An  H^^^'(kb)  +  H^2^’(kb)  =  icp  (14) 

since  the  cosine  functions  are  linearly  independent. 
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In  writing  the  general  solution  for  the  sound  pressure 
Pg  outside  the  shell  it  is  necessary  to  account  for  the  fact  that, 
for  large  r,  the  pressure  field  should  have  the  form  of  an  out¬ 
going  wave.  Since 


Lim 

r-'os 


H(l)(kr)|  =  Vs/t,  kr 


represents  an  outgoing  wave,  the  proper  form  for  p^  is 


p2  =  2  H^^^(kr)  cos  n0  .  (15) 
n=0 

The  interface  condition  on  the  velocities  of  the  fluid  and  shell 
on  the  outside  surface  of  the  shell  v;ill  be  matched  at  the  middle 
surface  of  the  shell,  as  in  Eq.  6.  Hence 


1  »Pa 


imo  d  r 


=  w 


or 


r=b 


^  Z  '  (kb)  cos  n9  =  cos  n8  , 


imp 

n=0  n=0 

and,  noting  the  linear  independence  of  the  cosine  functions, 


(16) 


Cn  H^’-^'(kb)  =  icp  .  (17) 

Equations  11,  14,  and  I7  give  three  sets  conditions  on 

the  four  sets  of  unknowns,  A„,  B  ,  C  ,  and  G  .  A  consideration 

n'  n'  n’  n 

of  the  motion  of  the  shell  is  necessary  to  obtain  the  requisite 
fourth  condition. 
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B .  Mo  tip  n  of  the  Shell 


Let  the  radial  and  tangential  displacements  from 
equilibrium  of  the  middle  surface  of  the  shell  be  denoted  by 

•7 

Uj.  and  Ug ,  respectively.  FlQgge*s  shell  equations  are' 


's  T 


_s 
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d  u. 


be‘ 


d0 


(18) 


Ps 


3u^ 

■iT 


u  ] 

r-' 


E  h^  3^u  dV 

- r  [ - F  +  2  — ^  +  u  ]  (19) 

I2b^  be  ae  ^ 


where  h  is  the  thickness  of  the  shell  and  b  is  the  radius  of  its 

middle  surface,  p.  is  the  mass  density  of  the  shell  material  and 
P  s 

E  =  E/ (1-v  ) ,  E  being  Young's  modulus  and  v  Poisson's  ratio. 

o 

p_  and  p^  are  the  pressures  on  the  inside  and  outside  shell 
surfaces,  respectively,  so  p_  =  P4.  =  [p2]r=b‘ 

Since  w  is  the  radial  velocity  of  the  shell. 


w 


=  -imu^  . 


(20) 


The  tangential  velocity  v  of  a  shell  element  is  represented  as 


V 


at 


=  -imu, 


(21) 


For  S3mimetry  about  6=0,  the  Fourier  series  representation  for  v 
is  of  the  form  ^ 


V 


^  Dn  Sin  n0 
n=0 


(22) 


The  shell  Equations  18  and  19  can  be  differentiated  with 
respect  to  time  to  give  conditions  on  u  and  v. 
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b2 


■UJ  pgW  = 


lU) 


[  Pi 


r=b 


]  -  -f  [ft 

r=b  b^ 


(23) 


[  A  +  2  ^  . 

ae  ae'^ 


(24) 


If  Equations  4,  l6  and  the  Fourier  series  for  v/  and  v  are 
substituted  into  Equations  23  and  24  and  the  linear  independence 
of  the  trigonometric  functions  employed,  the  following  conditions 
are  obtained. 

h  E  n 

-im  [H^^^Ckb)  A„  +  (kb)  B J +iu)  (kb)  C - §— 

ri  nri  n 

(25) 

(26) 

C .  Determination  of  the  Coefficients 

Equations  13,  15,  18,  25,  and  26  represent  a  set  of 
simultaneous  algebraic  equat 3.0ns  in  the  unknowns  A  ,  B  ,  C  ,  D  , 
and  G^.  This  set  of  equations  is  readily  solved  to  complete  the 
solutions  for  the  sound  pressure  field  in  Regions  1  and  2  given 
by  Equations  3  and  15  respectively. 

As  an  example  of  toe  form  of  the  coefficients,  the  solution 
for  the  set  of  coefficients  C.^  v?hich  are  associated  v/ith  x>2 
(Eq.  15)  is 


E  2 

+  h  p - i  Cl  + '’■^(n^-2n^+l)] 

'  ^  b  ia» 


n  G„  =  0 
n  ,2  n 
b 


G„  =  0 
n 
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c  =  , 

n 


n 


H^’(ka)  +  h  Hj;'(kb)  [J^(kb)  N^(ka)  -  N^^Ckb)  J^(ka)] 


(27) 


where 


7  = 
'n 


rrb 


2c% 


Ps  ■  % 

®  b  12b 


211*^  +  ?.)]  - 


r  2  „2 


-  Eg  n^) 


This  form  is  typical  of  the  remaining  coefficients.  It  is  of 
interest  to  note  that  for  h  =  0,  the  reduce  to  the  coefficients 
obtained  for  a  transducer  radiating  without  a  dome.  This  limiting 
case  can  be  compared  to  the  results  of  Appendix  where  the  bare 
transducer  is  studied  independently. 
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